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We calculate the first-order perturbation correction to the ground state energy and chemical 
potential of a harmonically trapped boson gas with contact interactions about the infinite repulsion 
Tonks- Girardeau limit. With c denoting the interaction strength, we find that for a large number 
of particles N the 1/c correction to the ground state energy increases as N 5 ^ 2 , in contrast to the 
unperturbed Tonks-Girardeau value that is proportional to TV 2 . We describe a thermodynamic 
scaling limit for the trapping frequency that yields an extensive ground state energy and reproduces 
the zero temperature thermodynamics obtained by a local density approximation. 
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I. INTRODUCTION 

The realization of quasi one-dimensional ultracold bo- 
son gases with tunable interaction parameters [1-3] and 
the succeeding advances in atom chip trap technology 
[4-7] have renewed interest in theoretical models of one- 
dimensional bosons with short-ranged interactions. Of 
particular relevance to workers in this field is the Lieb- 
Liniger model [8] in which contact interactions are de- 
scribed by Dirac delta functions. The suitability of this 
model in describing the low temperature properties of 
these quasi one-dimensional bosonic systems has been 
further strengthened by Olshanii's analysis of the low 
energy scattering of atoms under tight transverse har- 
monic confinement: The longitudinal s-wave scattering 
amplitudes are indeed reproduced by a one-dimensional 
pseudopotential proportional to a Dirac delta function 
[9]. The resulting demonstration that the magnitude of 
this effective delta interaction can be explicitly calculated 
from the three-dimensional atomic scattering length and 
the dimensions of the confining external trap [10] further 
strengthens the link between theoretical one-dimensional 
models and quasi one-dimensional experiments. 

Still, the free Lieb-Liniger model is quite an idealiza- 
tion for actual experiments because the atoms are gen- 
erally longitudinally confined by an external potential 
and thus much effort has been devoted to studying the 
effects of confinement of interacting bosons [9, 11, 12]. 
Introducing an external harmonic potential to the free 
Lieb-Linger model of spinless bosons leads to the many- 
particle Schrodinger eigenvalue equation 



h 2 d 2 



N 

^~2^i~dx 2 

,i=i l 



+ 



2 2 

mu! xf 



2 

ma 



]6( Xj - Xi )* h , (1) 



* Electronic address: fparaanOic.sunysb.edu 

t Electronic address: korepin0max2.physics.sunysb.edu 



where m is the mass of each of the N atoms, lo is 
the angular frequency of the trap, and a is the one- 
dimensional scattering length. The superscript 'b' refers 
to the bosonic nature of the labeled quantities. Measur- 
ing energy in units of hw and length in oscillator units 
£ = \J h/mui gives the dimensionless eigenvalue equation 
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where we have introduced the dimensionless interaction 
strength c = —Ija. We consider here the repulsive case 
c > (negative scattering length) to be specific. In 
the absence of an harmonic potential the correspond- 
ing eigenvalue equation is solvable by the Bethe ansatz 
and consequently much is known about the ground state 
and elementary excitations of this system [8] , as well as 
the properties of the various correlation functions at zero 
and finite temperatures [13]. However, for the important 
case of harmonic confinement an exact solution to this 
problem for general values of the interaction strength c 
is lacking. The sole exceptions are the two-particle case 
that is separable in relative and central coordinates [14], 
and the Tonks-Girardeau (TG) limit of infinite repulsion 
c — »■ +oo in which the exact TV-particle wavefunctions 
are absolute values of Slater determinants [15]. These 
results, especially the latter, are prototypical examples 
of the fermion-boson duality derived by Girardeau [16] 
and later generalized by Cheon and Shigehara [17] for 
one-dimensional systems of particle's having contact in- 
teractions. For finite values of the interaction parameter 
c, expressions for the atomic density and collective oscil- 
lation frequencies have been calculated using local den- 
sity approximations [9, 11, 18-20] and time-dependent 
density functional theory [20, 21], while formal expres- 
sions for the self-consistent Hartree-Fock equations for 
the single-particle density matrix have been obtained for 
general trapping potentials [12]. The analogous problem 
of a system of confined interacting fermions has also been 
treated by similar methods [22]. 

In this work, we will use the mentioned fermion-boson 
relation to develop perturbative 1/c corrections to the 
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ground state energy and chemical potential of an harmon- 
ically confined interacting boson gas about the Tonks- 
Girardeau solution. The details of the specific fermion- 
boson mapping we employ here that utilizes a non-local 
pseudopotential [12] are given in Section II. In Section 
III we present our perturbation calculations for the gen- 
eral case of an iV-particle system obtaining a closed form 
analytical result that is calculable for any N. We ana- 
lyze few-body cases and discuss the thermodynamic limit 
TV — > oo of our solution in Section IV. We summarize our 
results and give concluding remarks in Section V. 




II. FERMION-BOSON MAPPING 

In one dimension it has been demonstrated that a 
bosonic model with pairwise contact interactions of 
strength c can be mapped into a fermionic model with 
pairwise interactions of strength 1/c [17]. Specifically, 
given a fermionic wavefunction $ f that satisfies the eigen- 
value equation 
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an appropriate choice of a pseudopotential operator V { 
allows us to make the following correspondence between 
the bosonic wavefunction \E' b and its fermionic counter- 
part <I> f 



* b = A$ { (xi,...,x N ), 
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where A = Yli<j s S R ( x j ~ x i) 1S a function that is com- 
pletely antisymmetric under any transposition X4 <-> Xj. 
Since the complex squares of the wavefunctions ^f h and 
$ are identical, the energy eigenvalues E h and E are 
equal. For our specific problem, a suitable fermionic 
pseudopotential operator has matrix elements in the co- 
ordinate representation given by [12] 
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where 



= Xj — xi and Rij = \{xj + Xi) are rel- 
ative and central coordinates, and ty? f (cci, . . . , a: at) and 
4> f (xi, . . . , xn) are the coordinate space wavefunctions 
corresponding to the fermionic state kets \ip l ) and \cf> f ). 

In the infinite repulsion limit, the bosonic eigenvalue 
equation (2) is solved by the absolute value of the ground 
state Slater determinant [15] 
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where $> l TG is the fermionic ground state wavefunction, 
{x m } are the coordinates of the atoms and {ipn} are the 
N lowest energy single-particle harmonic oscillator eigen- 
functions ip n (x) = TT' 1 ^(2 n n\)~ 1 ^H n (x)e- x ^ 2 . The 



FIG. 1. The first-order perturbation result (dashed line) for 
the ground state energy of two delta interacting bosons in a 
harmonic trap is compared to the exact solution (solid line). 
The ground state energy -Etg in the infinite repulsion limit is 
given by the horizontal dotted line. 



H n (x) appearing here are the usual Hermite polynomi- 
als. The corresponding energy of this TG ground state 
(in units of hto) is 
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For a finite and large repulsion strength we may therefore 
use the quantity 1/c < 1 as a perturbation parameter 
for the fermionic problem (3) so that ordinary first-order 
perturbation theory gives the desired correction 
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III. GROUND STATE ENERGY CORRECTION 

In this section we work mainly in the fermionic sec- 
tor and omit the superscripts 'f for brevity. Here, our 
objective is to explicitly evaluate the leading correction 
SE = (^tgI^I'&tg) to the ground state energy Eq. Since 
$tg is a Slater determinant and the fermionic interac- 
tion operator V is a sum of two-body operators v having 
matrix elements 
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we may calculate the perturbation 5E using the Slater- 
Condon rule 5E = J2k<i( v kiki - vuik) [23]. Prior to 
calculating the derivatives appearing inside the integral 
(9), we must be careful to write the coordinates x\ = 
R+\r and x 2 = R—\r in terms of the relative and central 
coordinates r and R. The symmetry of the integrand 
allows us to write and define Vkm = —Vkiik = iu where 
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Thus, the leading correction becomes 5E = 2 X)fc<z ^kl 
and is always negative as expected. After some manip- 
ulation, we obtain a finite series expression for the en- 
ergy correction that may be evaluated for any number of 
atoms N: 
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IV. 



FEW PARTICLE RESULTS AND 
THERMODYNAMIC LIMIT 



The special case of N = 2 particles is separable in 
relative and central coordinates and the resulting eigen- 
value problem for the ground state energy is exactly solv- 
able. Upon imposing vanishing boundary conditions on 
the two-body wavefunction at infinity, we find that the 
ground state energy E of the trapped two boson system 
satisfies the transcendental equation 
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where Eq e [1,2]. Our pcrturbative result (11) gives 
5E = — (2/c)\J2/tt and coincides with the leading cor- 
rection term in the 1/c series expansion of this exact so- 
lution about the TG ground state energy E T q = 2. Also, 
an essentially identical result was obtained by Sen [24] for 
the reduced single-particle problem using a different lo- 
cal pseudopotential involving the second derivative of the 
delta function cx S"(x 2 — xi) [25]. A comparison of our 
first-order perturbation result and the exact two-particle 
ground state energy is shown in Figure 1 and we observe 
good agreement between the two results in the strongly 
interacting regime c 3> 10. 

Before discussing the situation for more than two parti- 
cles, let us restore units and rewrite the perturbed ground 
state energy as 

E (N) = \tuoN 2 [1 + 2a(N)/c] + 0{l/c 2 ), (13) 

where a(N) is a dimensionless function of N. For val- 
ues of N up to 10 3 the magnitude of the scaled first-order 
correction —c5E(N)/N 2 is plotted in Figure 2a as a func- 
tion of N on a double logarithmic plot. Inspection of this 
graph suggests a simple power law scaling for the first- 
order correction with large N. Leading order asymptotic 
analysis reveals that the partial sums in Eq. (11) scale as 
^AT3/2 for krge N so that TV 2 a (AT) ~ TV 5 / 2 : 

E (N) w \HujN 2 [l + 2a V A /V/c] , N -> 00, (14) 

where ao is a constant number. Indeed, for as few as N > 
15 particles the factor a(N) is quite well approximated 



(within 1%) by the function aov N with ao ~ —0.408. 
In other words, the correction factor a(N) reaches its 
asymptotic scaling behavior for systems as small as N ^> 

olio 1 ). 

To obtain a thermodynamic limit with an extensive 
ground state energy, we observe that in addition to send- 
ing the number of particles to infinity we must also re- 
quire the trapping frequency ui to vanish as 1/N. As 
we shall see, this condition allows us to reproduce the 
Thomas-Fermi results near the Tonks- Girardeau limit. 
This scaling requirement implies that the quantity £/\/N 
approaches a constant in the thermodynamic limit, in 
contrast to the thermodynamic limit in the case of a Lieb- 
Liniger gas confined in a flat-bottomed box, in which the 
linear dimension of the system is taken to scale propor- 
tionally with particle number [26]. Looking back at our 
asymptotic expression a(N) « ctoVN, we find that the 
quantity y/~N /c approaches a constant value as N — > 00 
in our prescribed thermodynamic limit. This is precisely 
the condition used by Ma and Yang [11] to obtain the 
ground state energy of the trapped interacting boson gas 
in the Thomas- Fermi formalism, which we reproduce here 
in Figure 2b. We find that our first-order 1/c result is 
reliable for \/N /c < 0.1, which means that for a typical 
experimental setup with hundreds of atoms first-order 
perturbation theory and the Thomas-Fermi result coin- 
cide only in the extreme repulsion limit c > 10 2 . 

If we now define the Tonks-Girardeau limit chemical 
potential as /xtg = li m w-s-oo,a;->o It^jN and the scaled 
interaction parameter as 7 = liniAr^oo.^^o c/y/N, we 
obtain the zero temperature chemical potential 
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The first term in this expression corresponds to the chem- 
ical potential of free fcrmions in a one-dimensional har- 
monic trap while the last term gives the reduction in the 
chemical potential due to the finite repulsion correction 
and serves as a measure of the departure of the system 
from the unitarity limit. 



V. CONCLUDING REMARKS 

In this work we have calculated the first-order finite 
repulsion correction to the ground state energy of har- 
monically trapped bosons having contact interactions for 
any finite number of particles N. For N ^> ©(lO 1 ) we 
found that for a fixed interaction strength this correction 
scales as a power law A 5 / 2 , which allowed us to describe 
a thermodynamic limit that reproduces known results 
from Thomas-Fermi approaches. This contribution clar- 
ifies the smooth transition of the ground state properties 
of an harmonically confined interacting boson system as 
the particle number goes from finite N to infinity near 
the Tonks-Girardeau limit. We have demonstrated that 
in this strongly interacting regime, to at least leading or- 
der in 1/c, finite number effects are negligible in current 
experimental situtations that have ~ 10 2 atoms. 
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FIG. 2. (a) The scaled perturbation —cSE/N 2 (dots) grows as a power law y/N in the limit N — > oo. (b) Our large N result 
is consistent with numerical calculations (solid line) in the Thomas- Fermi approximation for y/N /c«l. The dashed lines are 
regression fits c5E ~ — 0.4087V 5//2 calculated from values N g [100, 1000]. The prefactor —0.408 is independent of the particular 
values chosen in this interval. 



A natural extension of this work would involve higher 
order corrections to the ground state energy and many- 
body wavefunction, as was done recently for a wedge- 
shaped trapping potential [27]. If we take the set of 
all fermionic Slater determinants as an expansion basis 
for ordinary perturbation theory about the TG limit, we 
discover that the perturbing pseudopotential couples the 
ground state to an infinite number of excited states. We 
therefore expect a complicated analytical result for the 
second order energy correction resulting in a numerical 
problem that may require a truncation of the correspond- 
ing Hilbert space. However, on the basis of the agreement 
between our asymptotic results and the Thomas-Fermi 



calculation (Figure 2b), we conjecture that the second 
order correction scales as N 3 /c 2 > in the thermody- 
namic limit. 
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